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The intensity of classical bright solitons propagating in linearly coupled identical fibers can be 
distributed either in a stable symmetric state at strong coupling or in a stable asymmetric state if 
the coupling is small enough. In the first case, if the initial state is not the equilibrium state, the 
intensity may switch periodically from fiber to fiber, while in the second case the a-symmetrical 
state remains forever, with most of its energy in either fiber. The latter situation makes a state of 
propagation with two exactly reciprocal realizations. In the quantum case, such a situation does not 
exist as an eigenstate because of the quantum tunneling between the two fibers. Such a tunneling 
is a purely quantum phenomenon which does not not exist in the classical theory. We estimate the 
rate of tunneling by quantizing a simplified dynamics derived from the original Lagrangian equations 
with test functions. This tunneling could be within reach of the experiments, particularly if the 
quantum coherence of the soliton can be maintained over a sufficient amount of time. 

Lead Paragraph 

Usually solitons in optical fibers are assumed to be classical (= non quantum) objects 
because they are made of a large number of photons. Nevertheless there exist quantum 
effects without classical counterpart, like the tunneling under a potential barrier. We 
investigate one possible realization of such a quantum tunneling with solitons as basic 
entities. Specifically, we consider a soliton propagating in two linearly coupled fibers 
that are assumed identical. It has been known for some time that, at small enough 
coupling, asymmetric solitons only can propagate and be stable. The amplitude of 
such asymmetric solitons is predominantly in either fiber and remains there forever 
classically. This makes, for a given energy, two possible steady states exactly symmet- 
rical with respect to each other under permutation of the two fibers. In the quantum 
version of the same problem, the two solitons merge into a single quantum state shar- 
ing a quantum amplitude spread between the two fibers, because of the possibility of 
quantum tunneling from one fiber to the other. We study this problem thanks to a 
reduced set of equations derived from the full set of coupled nonlinear PDE's by choos- 
ing convenient trial functions for the classical soliton dynamics. Thanks to this choice, 
the bifurcation pattern of the soliton solution in the coupled fibers is well reproduced. 
Because the trial set generates dynamical equations with a Lagrange structure, this 
Lagrangian system is relatively easy to quantize. To obtain the quantum amplitude of 
transmission by tunneling under the barrier, one replaces the original Hamiltonian sys- 
tem by its Euclidean counterpart. Orders of magnitude relevant for a possible physical 
application are given. 

PACS numbers: 

I. INTRODUCTION 

Generally speaking a soliton is a localized solution of an equation for the propagation of a field envelope. It stays 
localized under the opposite effects of linear dispersion tending to spread the wave and of nonlinearity making the 
wave steeper. We make one more step by considering this soliton as a 'true' particle, that is by seeing it as a classical 
object that should be ultimately quantized to keep the consistency of our view of the physical world. This is of course 
not a new idea, see for instance the review [l[ on the quantization of various nonlinear equations for classical fields. 
Quantum effects are irrelevant for macroscopic phenomena like solitary waves in a water channel. However, there 
is an instance of solitonic physics where quantization could bring significant new effects, namely the propagation 
of optical solitons in fibers: there the amplitude of the wave may be small enough to yield solitons with not too 
large action, measured in units of Planck's constant Jr-. If the action is much larger than this quantum unit one is 
in the classical regime, many trajectories contribute to the saddle point of the Feynman integral [2j and quantum 
interferences between coherent quantum states become practically impossible. Conversely, if the action is not too large 
compared to this quantum unit, one could observe quantum phenomena as tunneling and interferences. Moreover, 



even if the quantum state has a coherence time shorter than the tunneling time, there is still quantum tunneling, 
but at a reduced rate because the build-up of the state on the other side of the barrier is slowed down . Below 
we show that quantum tunneling of a soliton may occur between two weakly coupled fibers, and we discuss the 
possibility of quantum interference between the two states carried by each fiber. The starting point of our study is 
the well established result that classically a soliton injected in a given fiber cannot switch to the other fiber, when 
the coupling is less than a certain critical value. In this range the soliton evolves towards the stable a-symmetric 
solution having its energy predominantly in the initial fiber [i| . There exist another stable solution which is obtained 
by permutation of the two fibers. At small coupling the two a-symmetrical states are separated by a finite barrier 
that cannot be crossed classically. We predict that this may be wrong in practice because of quantum tunneling. 
Our derivation is based on the calculation of the tunneling probability which writes in the WKB approximation, as 
T = exp(—2S/h), where S is the physical action associated to tunneling [5(. The quantity S/H may be small, even 
for a pulse with a large number of photons, because all dynamical phenomena we consider, like the balance between 
the nonlinearity, the group velocity dispersion and the coupling, imply small perturbations to the dominant effect 
resulting from the linear dispersionless terms of Maxwell's equations. The perturbations we consider are as small as 
the nonlinear term 122! with respect to the dominant term no, in the expansion of the refractive index n = hq + n^l 
for the Kerr medium of the fiber, / being the optical intensity. In this sense a soliton is a bound state of photons: 
photons are attracted to each other by the focusing nonlinearity. The soliton resembles the atom of a heavy element 
which is a quantum object made of many electrons and nucleons. Atomic physics has also to do with energies much 
smaller than the rest energy of the particles (electrons and nucleons) making the atom. Another idea of atomic (and 
quantum!) physics is relevant for our goal: If the atom remains in its ground state, it is able to make interferences 
with a wavenumber depending on its mass and velocity only, independently of the details of the state of its electrons 
and nucleons. In our study we assume the absorption and change of frequency of the photons by inelastic and/or 
Raman scattering to be negligible, and discuss the role of these effects in interference experiments in the last section. 

Let us outline the organization of this paper. In the Section |TT] we introduce the classical model of propagation of 
solitons in fibers. First we deal with the single fiber, then with the two coupled fibers. There is nothing new here and 
we focus on what is relevant for us, namely the bifurcation diagram as a function of the coupling. For strong coupling 
the stable solution is symmetric (with intensity equally shared between the two fibers). Below a critical coupling, the 
classical prediction is that the stable soliton is a-symmetric, with a high amplitude in one fiber and a small amplitude 
in the other, as said above. The details of the transition are a bit complex because it is subcritical. This has been 
studied [H by direct numerical solution of the coupled PDE's describing this problem (equations ((U and ((TJ below). 
However interesting it is, this model presents some difficulties for our quantization problem. Therefore, in the next 
section IIIII we outline another approach to the same problem, namely we use the Lagrange formalism to compute 
approximate solutions with trial functions (instead of the full unknown solution) . Thanks to an approriate choice of 
these functions, the pattern of bifurcations of the asymmetric to symmetric solitons, known from the direct numerical 
simulations, is well recovered. We use the same trial functions as Malomed et al. [y] who discuss very thoroughly the 
general issue in a paper that we recommend to the interested reader. 

As explained in section [Tvl this 'trial dynamics' is used then to quantize the system. This is of course not exact, 
but requires far less formalism than the full quantization of the two coupled nonlinear field equations. Thanks to this 
method, one can use standard results and methods of quantum mechanics for systems with a few degrees of freedom 
(as opposed to field theories). In particular, we can compute the trajectory under the potential barrier, found by 
multiplying the propagation variable z, and the Hamilton- Jacobi action S by i. This yields a well defined problem 
of Hamiltonian mechanics, called sometimes the Euclidean version of the initial problem. The tunneling factor is 
derived from the action of the heteroclinic trajectory joining the two equilibria: a stable equilibrium in the original 
Hamiltonian system remains an equilibrium in the Euclidean one, but it becomes unstable there. It turns out that 
the tunneling probability depends algebraically on the coupling between the two fibers. This is a significant remark 
for possible applications, because it yields a much smoother dependence with respect to the coupling than the usual 
exponentially small tunneling amplitudes. 

The last section summarizes the main results of this paper, discusses the possibility of interferences and presents 
some ideas for possible applications. Quantitative predictions rest on the rather complex problem of turning back 
from the dimensionless equations used throughout this work to quantities with a physical dimension, something done 
in the Appendix. 



II. CLASSICAL PROPAGATION OF SOLITONS IN COUPLED FIBERS: THE GENERAL MODEL 



The mathematical model for the dissipationless propagation of optical solitons in one fiber is the classical (= non 
quantum) nonlinear Schrodinger equation: 

dE d 2 E , , 2 
l ^ +a ^ +m E = °- (1) 

Even though this equation is called nonlinear Schrodinger (NLS), it does not mean at all that it makes a quantum 
system. It only resembles the usual Schrodinger equation, but it describes a purely classical field, exactly as Maxwell's 
equations do for an EM field. This equation is written with real coefficients a and (3 carrying a physical dimension 
to make possible the discussion (see Appendix) of the order of magnitude of the physical effects to be expected. The 
field E is the complex amplitude of the electric field in the wave. We take it as a scalar, although polarization effects 
could be brought into the picture in principle. This equation is derived in the Fresnel approximation, assuming that 
the changes of amplitude along the fiber are much slower and on much longer scales than the oscillations of the optical 
field itself and it is also written in the frame of reference moving with the speed of the envelop of the wave, where 
the position variable is z. For a and (3 real, this equation has a Lagrange-like structure. It cancels the first order 
variation of the 'action' 
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In this equation E is the complex conjugate of E. The writing of the action in equation ^ brings in an important 
problem, because it is not 'the' physical action. Such a physical action has to have the dimension of the product of 
an energy and of a time. Therefore, the action written in equation cannot be an action from the point of view 
of physical dimensions. The physical action of the EM field is proportional to S, its derivation is postponed to the 
Appendix. An overall constant multiplying factor does not change the Euler-Lagrange equations, but it is crucial 
when quantizing this system because this relies on a comparison between the action and h, two quantities with the 
same physical dimension. 

By rescaling E — > ii/3 -1 / 2 (assuming (3 positive to be in the focusing case where solitons exist), and t — > t(2a) 1 / 2 , 
one obtains the dimensionless nonlinear Shrodinger equation: 

.du 1 d 2 u . ,, 

*dz + 2W + M U = °- (3) 

This equation has a number of interesting symmetries. In addition to the Galilean invariance (if u(z, t) is a solution, 
then u(z, t — 2;/C)e^^ t_ is also a solution), it has a dilation symmmetry: if u(z, t) is a solution and \x an arbitrary 
real number, then /xit(z/it 2 , tjj,) is also a solution. It has a two parameters family of soliton solutions: 

ue i(p 

u = — — r = ve ltp sec\i(vi). (4) 

cosh(i/i) 

In this solution, v is any real number and the phase ip is if = \v 2 z + ipo, with ipo arbitrary constant phase. 

Note that z plays here the same role as the time in the usual Schrodinger equation. Among the conserved quantities 
associated to any solution of the NLS equation, let us write the " energy" 
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Suppose now that, instead of a single optical fiber, we have two identical coupled fibers, and that the coupling is 
linear and preserves the symmetry between the fibers. The propagation of solitons in this system has been studied in 
the last fifteen years 0], 

To describe the two coupled fibers supporting solitons we introduce two focusing NLS equations, written in a 
dimensionless form: 



du 1 d 2 u 2 



-KV, 



and 



.dv 1 d 2 v 



dz 



2 dt 2 



+ \v\ 2 V — —KU, 



(6) 



(7) 



where k is the strength of the linear coupling, and we define the "mass" 

Q = J(\u\ 2 + \v\ 2 )dt. (8) 

which is a constant of motion. Consider solutions of the form u(z,t) = U(t)e lqz and v(z,t) — V(t)e lqz . Because of 
the common phase factor e lqz the ^-dependence cancels out and the two functions U(t) and V(t) are solutions of the 
two coupled ordinary differential equations: 



-qU+^ + \U\ 2 U=- K V 

-qV + ^ + \v\ 2 v = - K u. {> 

For the solution to decrease to zero when t tends to plus or minus infinity one must have q > 0. Furthermore the 
sign of n can be changed by changing U into —U for instance and keeping V the same. We choose k positive that 
corresponds to in-phase stationary solutions (U, V), the out of phase ones being unstable [4j. 

This set of equations has been studied numerically and analytically Q. An exact calculation shows that the 
symmetric solution U = V always exists and is linearly stable in the range -S= < -7=. For higher values of this 
ratio, the symmetric solution looses its stability and an a-symmetric solution branches off. While it is not explicitely 
mentioned in [H, the subcritical character of the bifurcation can be deduced from Fig. 11 of the paper by Akhmediev 
and Soto-Crespo(1994) when using as control parameter and ^ as order parameter, i.e. by rotating the figure. 

Consequently no stable and weakly asymmetric solutions branches off the unstable symmetric soliton for k slighty 
smaller than the onset of linear stability, although an unstable asymmetric solution branches off at values of k slightly 
larger than the critical one. Furthermore a branch of stable asymmetric solitons goes continuously from k = to a 
finite coupling, slightly larger than the value of linear instability of the symmetric soliton. The stable asymmetric 
soliton disappears by a saddle-node bifurcation for a value of the coupling that is, by a numerical coincidence, very 
close to but smaller than the onset of linear stability of the symmetric soliton. At this saddle-node bifurcation the 
unstable and stable a-symmetric solutions merge to disappear at smaller values of 

In the numerical investigations of this problem an interesting phenomenon comes into play, namely the radiation of 
energy at large distances of the solitons. The amount of radiation is stronger when the initial conditions are further 
away from a stable solution [4j. This radiation happens in the far wings of the time dependent amplitude profiles 
(|ti(£)| , |w(i)|), where the full equation reduces to its linear part. Although very strongly dispersive this describes 
radiation by wave packets of ever increasing width, but carrying nevertheless energy and eventually other invariants 
to infinity. Such a coupling between a localized system and the infinitely many degrees of freedom of a radiating 
field may lead to irreversible effects [ll|. It shows how subtle may be the distinction between 'dissipative' and 
'nondissipative' systems as soon as one goes beyond the obvious. Irreversible process due to radiation may not even 
require an infinitely extended physical space. They may also take place in the reciprocal (or momentum) space by 
cascade of energy toward smaller and smaller scales, a typically nonlinear phenomenon We plan to come back to 
the issue of the effect of radiation on quantum phenomena in the present problem. We shall neglect this kind of effect 
in the following, since they cannot be taken into account within our simple formalism. Even though the radiative 
losses are present, it was shown by Fadeev and Korepin [1] that they do not destroy the solitons in a single fiber, 
when they are included in the quantized version of the NLS equation. 

In the coming section we shall derive a reduced set of equations describing the propagation of soliton in coupled 
fibers. Indeed this reduction from the original PDE's to a set of coupled ODE's cannot be quantitatively exact. 
However with the same choice of trial functions as Malomed et al. G] we obtain at least a reduced system with the 
right pattern of bifurcation at decreasing coupling. The fundamental interest of this reduction is that it allows us to 
quantize the dynamical system rather straightforwardly. 



III. CLASSICAL PROPAGATION OF SOLITONS IN COUPLED FIBERS: THE REDUCED DYNAMICS 



Because of the lack of analytical solution in general, we follow an idea used already by various authors, that allows 
to understand in a fairly detailed way the results of the direct numerical simulation by using an analytical approach. 
This follows the general method of research of extrema of functionals by trial functions: dynamics can be reduced 
to a minimization problem, then one restricts the function space where this minimization is done to a space of trial 
functions depending explicitly on a few parameters and one studies the dynamical properties in this reduced space. 
Since we know the results of the direct numerical simulations it is in principle possible to check the quality of the 
approximation by comparing its predictions and the 'exact' results. This is necessary because the method of trial 



functions does not rely on a small or large parameter and so cannot hope to be 'exact' or close to exact in the usual 
mathematical meaning of the word. The papers by Malomed et al. @ discuss in depth the choice of the trial functions. 
We shall not reproduce this discussion here where we take their set of 'optimized' trial functions, following as much 
as possible their notations. 

The starting point is the writing of the action for the coupled NLS equations: 
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As it can be checked the action Snls is proportional to (z 2 — Zi) whenever the functions u{z,t) and v(z,t) are 
stationary solutions (with respect to the variable z) of the two coupled NLS equations j6| and (0 or functions 
proportional to the same phase factor e lqz . The problem we consider now is how does the coupling change the 
propagation of solitons. For that purpose we reduce the dependence with respect to t to an imposed form (the trial 
function) with arbitrary z-dependent coefficients, the trial parameters. Putting this trial form into the action integral 
and performing the integration over the variable t yields a functional of the parameters of the trial function that are 
themselves functions of z. Doing now the variation with respect to those functions, one finds at the end a set of 
ODE's for functions of z only. 

The choice of the trial functions is inspired by the soliton solution in a single fiber and it respects the symmetry 
between the two fibers. Following Uzunov et al.[6j one takes: 



u(z,t) = a{z) y/r](z)sech [rj(z)t] cos(Q(z))exp [i ($(z) + #(z) + q(z)t 2 )] , (11) 

and 



v(z, t) = a(z)y/^{z)sech [r)(z)t] sin(6(z))exp [i ($(z) - *(z) + q(z)t 2 )] . (12) 

In the case of a single fiber carrying a soliton of amplitude u, the trial function u(z, t) becomes the exact one-soliton 
solution with a = ^/rj, = q = 0, $ = zrj 2 /2 and constant. It is important to notice here that the angle O is for 
describing the balance between the two fibers, although the angles $ and \& have a physical meaning independent on 
the trial function, being related to the phase of the functions u and v. The angle could be replaced by another 
parameter in another trial function, not necessarily a circular function. Inserting this trial form into the action Snls 
and performing the integration over t, which is possible because the dependence with respect to t is fully explicit in 
the trial functions, one finds a reduced action that is itself the integral over z of the Lagrange function: 



C = 2^(z)cos(2*)sin(2e)-2a^cos(2e)- -a 4 r? sin^(29) + -a 4 ?? - -afrf - 2a z - — [-± + 2q*]. (13) 

dz 3 3 3 dz 6rj \dz J 

Up to obvious change in notations (from our z to £, from n to K, etc.) this Lagrange function is identical to the one 
written by Uzunov et al. [6( but for a misprint in their paper where the term q 2 in the last parentheses became g 4 
without harming the rest of their calculation. The parameters of the trial function are five functions of z: a, O, 'J', 
rj and q. The equations of motion for those five functions are derived by variation of the action, namely the integral 
over z of C. They read: 

derived by variation with respect to and 



' 9jf = -2Ksin(2*) 

sin(26)ff = ^77sin(26)cos(2e) - kcos(2*)cos(26) , , 

H - - V + £ [V 4 - aW (1 - isin 2 (29))] 

derived by variation with respect to 'J, 0, q, and rj respectively 

The parameter a 2 can be absorbed in the redefinition of k and will be set to 1 below (that corresponds to a mass 
Q 2i. 

The soliton solutions are z-independent solutions of this set of equations. There are two classes of soliton solutions 
in this model, depending on the coupling parameter. For any coupling there exists a symmetric soliton, with equal 
intensity in both fibers, i.e. O = At small coupling this symmetric solution is unstable against asymmetric soliton. 
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FIG. 1: Inverse of coupling coefficient k 1 versus the stationary values of x — cos(2Q). For large coupling, i.e. k > 1/6, the 
symmetric solution, x = 0, is the only one stable. 



Such an asymmetric soliton is found by canceling the z-derivatives in equations (|15p and choosing cos(2 l F) = 1. This 
yields the relation between the coupling coefficient k and the balance parameter for the intensity x — cos 20 

k- 1 = 6 - (16) 

which is illustrated in Fig.l. Using the trial functions, the bifurcation in the set of possible solutions is found to 
occur at the critical value, k c — 1/6 which is 11 per cent less than the exact value K c = V^3/4 . Moreover the 
bifurcation is slightly subcritical, in good agreement with the NLS results [||. Close to the bifurcation point there 
is a small range of values of the coupling, ^ < k < (0.167 < k < 0.181) where there are three sets of solutions: 
the symmetric solution that is linearly stable, and two pairs of asymmetric solutions, one linearly stable and another 
linearly unstable. The branch of stable asymmetric solution does not merge smoothly with the symmetric solution, 
but disappear when it has still a finite amplitude. The main conclusion that we shall draw here is that this set of trial 
functions reproduces well the pattern of bifurcation of the exact model. This makes it a good candidate for studying 
the quantum tunneling. 

Before to start this study, let us explain how we managed to define a quantity related to the usual potential energy of 
a mechanical system. Although this is not strictly necessary it helps to draw various quantities relevant for analyzing 
the tunneling by making a connection, however loose it is, with the familiar notions of barrier and of barrier crossing. 

The 'potential energy' is derived from the total energy associated to the dynamical system under consideration, 
namely the equations (fTi|) to (fT"5)) . An expression for this energy is given by Uzunov et al. With a 2 = 1 it becomes: 

< al = -2 K cos(2vI/)sin(2e) + i„ sin 2 (26) - ^ + ^ + (17) 

Note that the Lagrangian (fl3|) includes terms linear with respect to first derivatives (with respect to 2). It means 
that the two successive operations of choosing trial functions and averaging over the retarded time t, lead from the 
Lagrangian formalism to the Hamiltonian one, with 
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The first term in the r.h.s. of equation (fT5|) will be the one responsible for the Euclidian action derived in the next 
section. As already noticed by Uzunov et al., equations (fl4|) to (fT5]l are the Hamilton equations of a two-degrees of 
freedom system. The two pairs of conjugate variables are ,2a; = 2cos(20)} , and {q, y = , i.e. the phase 

and amplitude differences, as well as the chirp and width, respectively. We are interested in the value of H for 
steady states, that turns out to be a simple function of the coupling k. That should give an idea of how the energy 
changes when the variables are different of their values in the steady state(s). In order to preserve the connection 
with a potential energy in the usual sense we impose that, at the equilibrium points, this 'potential' energy is at an 
extremum. This is realized (probably not uniquely) by plugging into H. the values of q and 'J at the various equilibria, 
that is ^ = and q — to cancel the conjugate momenta. This yields: 

Hpot = -2Ksin(29) + ^r; sin 2 (26) - ^ + \-. (19) 



This 'potential' energy depends on two parameters, x = cos(26) and 77, and it is plotted in Fig. 2 for various 
coupling strength to show the bifurcation of the equilibria from a single equilibrium at large coupling, Fig. 2(a), to 




FIG. 2: Level lines of the potential in the plane (x,rf), (a) for k = 0.5, only the symmetric solution x — is stable, (b)for 
K — 0.18 which belongs to the subcritical domain in Fig.l where symmetric and a-symmetric solutions are stable, (c) for « = 0.1, 
only two asymmetric solutions are stable 

a more complex pattern, as the coupling decreases. In particular, some sort of barrier is evident in Fig. 2(c). It 
separates the two deep minima of the potential lying each in the vicinity of (±1,1). Each minimum corresponds to 
one of the stable asymmetric soliton, although the unstable symmetric soliton at (0,1/2) is a saddle point of the 
potential energy. 

The above picture illustrates the known results: classically there is no way for a soliton initially in a given fiber, to 
escape through the other fiber, at low coupling, because of the barrier. Before to present the quantum version of this 
problem, let us precise what is the low coupling range in terms of physical quantities. Note first that the low-coupling 
range writes 

— > 6, (20) 

K 

for an incident soliton of the form u(0,t) = a/ch(at) injected in one of the two fibers (a = 1 above). Secondly 
let us define the scaling quantities in equations©-®, by using the soliton units, z = z phys /Lp, t = t phys /r , 

u, v = \J 2 ™ 2 L~dBi 2, and k = K phys Ln , where Lu = is the dispersion length, and B\ t 2 are the slowly varying 
amplitudes of the electric field (see appendix). The relation PU|) becomes 

n 2 I M =n' 2 \B 2 M \ > -Ao/t^, (21) 

7T 

or, u 2 Im > fx when introducing the switching length L c — 2K p hys defined for the CW linear regime. Using the 
relation (|A3|) . the low coupling range also writes 

L c > 3irL D . (22) 

IV. SEMICLASSICAL QUANTIZATION OF THE COUPLED FIBER SYSTEM 

Before computing the quantum tunneling, let us recall the main differences between the classical solitonic solution 
and its quantized form. In quantum mechanics a state localized on one side or on the other only is not an eigenstate 
of the system, because of the possibility of tunneling. Therefore if one starts at 'time' zero with all the amplitude on 
one side (meaning all the probability in one of the two possible asymmetric states), after the time of tunneling this 
will be transferred to the other side and eventually oscillate between the two sides. It is also possible to inject at the 
input of the dual core fiber, the quantum ground state, which is symmetrical, with equal amplitude in the two sides. 
To have a physical image of the process by which the transition occurs between the two asymmetric states, one may 
recall that the number of photons is not fixed in the quantized soliton, so that it fluctuates in both fibers. Therefore, 
the fluctuations may bring one fiber into the soliton state, although the other goes to the state without soliton, and 
the two states switch in the course of time, as studied below. 



Let us now outline how to compute the quantum tunneling between the two fibers. Because we have a classical field, 
the quantization of the coupled equations §E§ , © for the two fibers belongs to the general problem of quantization 
of field theories. Although this may be done formally, it requires a rather heavy machinery in any case. Fortunately 
there are various possible short cuts in this derivation. The most obvious one is to reduce PDE's system to a set of 
ODE's, by using trial functions depending on a certain set of unknown parameters. By refining the choice of trial 
functions ad infinitum, namely by introducing trial function with more and more parameters, one should converge in 
principle toward the exact result. But we will merely use the above described trial functions. The Euler-Lagrange 
condition of stationarity of the action yields a set of dynamical (in 'time' z) equation, that can be formally quantized 
because it has a symplectic structure. 

This is what we are going to do, except for one point. It is possible to short cut all this explicit quantization in the 
WKB limit, where the wave function is expressed by means of the classical Hamilton- Jacobi action, $ = Aexp(iS/h). 
This is the well-known quasi-classical limit, that restricts oneself to situations where any action involved is typically 
much bigger than h . This seems a reasonable limit, but it does not necessarily cover all possible situations- we shall 
come at the end to what seems to be 'the' standard experimental situation in this respect. The WKB limit is especially 
convenient for treating tunneling problems, because it amounts to calculate the imaginary part of the action (which 
is complex) and to put at the end h at the right place. Indeed the tunneling factor is is given by T — ea:f>(— 2/Se/?j), 
at leading order. Here Se is the imaginary part of the action, which enters then in the modulus of the wave function 
as a real exponent (instead of the usual imaginary exponent relevant for the classical limit of quantum mechanics). 
This imaginary part of the action is calculated by two steps. First one has to change the conjugate variables (q,p) 
into (q,ip) in the classical Hamilton-Jacobi formulation of quantum mechanics, the Hamiltonian H(q,p) becoming 
H{q, ip). Secondly one is left with a problem of extremalization of a new action, the Euclidean action, that is formally 
another problem of classical mechanics. For instance in the often presented problem of a particle of energy E in a 

2 

double well potential V(q), with Hamiltonian H = + V(q), the Euclidean action is calculated with the abbreviated 
action (l3j 



S E = / Pdq (23) 

Jq[0] 

2 

derived from the Hamiltonian He = —-§^ + V(q) but with the same energy as the one of the classical motion. For 
the potential this means that it gets rotated by 180 degrees, thus exhibiting two "hills" of maximal energy. The values 
of q[0] and q[zf] in equation ([23"]) are those of the classical turning points defined by E — V(q) . To calculate Se, one 
has to find a trajectory joining these points, namely to calculate an Euclidean path integral. This is performed by 
solving the Hamilton equations for the Euclidean Hamiltonian 
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(24) 



by taking as initial conditions, the known value q[0], and an unknown value p[0]. By varying the latter value, one 
finally converges towards a trajectory ending at q[zf], which provides the action defined in equation ([25]) . Note that 
equations (|2"4"|) are obtained from the classical Hamiltonian system ( which is identical to equation (|24[) but with H in 
place of He), by changing z in iz, and p in ip. The change to an imaginary "time" ( from z to iz here) amounts to 
go from a Minkowskian to an Euclidean metric. Therefore equations (|24|) are called "Euclidean equations of motion" , 
and their classical solution joining the two "vacua" of the double- well potential, often named "kink solution", is an 
example of an instanton [14j in quantum mechanics. 

In the above example the variables (p, q) are the impulsion and position of a particle in a ID potential. General- 
ization to cases of a multidimensional set of generalized coordinates and momenta leads to similar relations [l3[ . 



A. Semi-classical Action. 



To put all those principles in practice we have to formalize the dynamical system (equation (|15|) ) in terms of 
canonically conjugate variables. Once this is done, the Euclidean equations of motion are found by multiplying the 
"time" z and the momenta by i. As noted in section III, the reduced equations (jT5J) are those of an Hamiltonian 
system with two degrees of freedom, therefore a simple choice for conjugate variables (qj,Pj) with j — 1, 2, is to take 



the pair (x — cos(20),y — g^y) as coordinates and (2\E f ,g) as their conjugate momenta. The Euclidean Hamiltonian 
is obtained from the classical one in equation p9|) . by changing cos(2^) into cosh(2\E'), and q 2 into — q 2 . It becomes 

WE.trial = -2k cosh(2*) sin(26) + ^77 sin 2 (26) - ^ + ^ - (25) 

The semi-classical dynamics is then driven by the new set of four (Euclidean) equations, that are the Hamilton 
equations for the conjugate variables (qj,pj), deduced from the Euclidean Hamiltonian (J35j) 



' f£ = -2Ksinh(2*)Vl-2; 2 
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FIG. 3: Quantum trajectory superposed to the potential for k = 0.1: only two asymmetric solutions are stable. 

As in the case of a particle in a double well potential, calculating the probability for the soliton to tunnel through 
a classically forbidden region ( H pot ) with the Minskowskian space path integral, corresponds to calculating the 
transition probability to tunnel through a classically allowed region ( —H pot ) in the Euclidean path integral, with 
the action 



St 



(pidqi + p 2 dq 2 



(27) 



where qx, qn are the coordinates of the turning points. To perform the integration giving the action, it is enough to 
choose a convenient integration path in the Euclidean plane connecting the two minima M(xm ,Vm) of the classical 
potential in Fig. 3, which become the maxima of the Euclidean potential. For small values of k, one has 7i po t,M ~ 
— ^(1 + 18k 2 ). In the present case it is easier to carry the integral from xq — up to xm- We set the value of the 
Hamiltonian TIq close to 7i po t,M- Because of the symmetry of the heteroclinic trajectory joining the two extrema of 
the potential, we have to choose the initial condition qo = 0. Then the initial value of the phase difference ^0 is 
deduced from equation (|25[) . and we have only one initial parameter to adjust, 770, in order that the trajectory ends 
with a vanishing impulse qf = tyf — close to the extrema M(xm,Vm) in the plane (x,rj). The integration path is 
shown in Fig. 3. The action along the semi-classical trajectory is given by the expression (|27|) that writes with our 
notations 



Se = 2 



z=0 



-4*(z)Ksinh(2*(z)) v /T — x(z) 2 — 



2n 2 q 2 (z) 



3 rf{z) 



(28) 



The numerical result of the integration is shown in Fig. 4 which displays the action as a function of the coupling 
parameter k in a logarithmic scale. In the domain of existence of the asymmetric solution, k^ 1 > 6 in Fig. 1, the 
action clearly displays a logarithmic dependence with respect to the coupling, we have the law 



Log, (E 2 /K) 



FIG. 4: Action/2 (squares), Maximum of the impulse phase $o (circles), and 770 (triangles) versus the scaled coupling parameter 



S e = 2]h(k c /k) (29) 

that holds true with a precision better than 1 per cent over many decades, with the numerical value k c = 0.2 
slightly higher than the bifurcation one 0.18. Note that the relation (|2"9"|) holds true except in the close vicinity of 
the bifurcation point, not visible in Fig. 4. The Zn-dependence in equation (|29|) follows straightforwardly from the 
substitution of exponentials for the hyperbolic sine in the equation of motion for Euclidean dynamics. We also report 
in Fig. 4 the dependence of 'J'o and ?7o a s function of k. At x — the solution becomes transiently symmetric, 
sin(26) = 0, but its width is different from the symmetric value, 770 7^ |, and the impulse ^ is maximum. We show 
that \I/o evolves very much as the action, while 770 is quite constant. Actually the heteroclinic trajectory drawn in 
Fig. (3) passes through the abscissa x — approximately at the ordinate r\ ~ 0.67 whatever the value of the coupling 
constant, while the impulse here increases like ln(fc c /«;). This result shows the leading role of the conjugate variables 
*S> and x — cos(20) in the dynamics. At this stage it is interesting to compare the latter result f2"9"|) with the action 
derived by a simpler choice of trial functions, based on the hypothesis of constant width soliton (and of no chirp), 
as proposed by Pare Q and Kivshar (oj. In these simpler cases, one obtains a single degree of freedom Hamiltonian 
dynamics. The approximate calculation of the Euclidean action may be done analytically, and leads to similar results 
in both cases. With the notations of Kivshar, for example, using as conjugate variables ($, A), the calculation of the 
action amounts to carry the integral Se = $(A)dA, the function $(A) being given explicitly in 0. In the limit 
of a small coupling and with k — 7 -1 , the equation for $ reduces, at leading order, to <& « i 111(7) + ^ where <3E> is the 
solution of j^-y = e 4 * that is of order 1. Therefore in this limit 7 large (equivalent to small coupling), <f> w iln( 7 ) so 
that the action associated to tunneling is just S ~ 2iln(^-), where j c is a constant. 

Summarizing the Euclidean action obeys the law (|29|) in all cases of trial functions we have considered, i.e. for a 
single degree of freedom Hamiltonian as well as with two degrees of freedom. Consequently, it does not seem necessary 
to refine more our model to obtain the information we need, i.e. the order of magnitude of the tunneling amplitude. 

B. Tunneling factor. 

The possibility for the soliton to tunnel from one fiber to the other in real space, is measured by the transmission 
coefficient, with the expression T = | ^ | in a double well tunneling problem, with F : A the amplitudes of the 
transmitted and incident waves, respectively Q . It has already been noted that the transmission is given by T = 
exp(—2S/H) at leading order. In practice S is the "physical action", having the same dimension as H. Therefore, to 
calculate the "true" transmission for the soliton in the two coupled fibers one has to multiply the dimensionless action 
Se by an appropriate coefficient depending on the properties of the fiber and of the characteristics of the EM 
wave, this giving lastly the "physical action" S E hya — s^Se which has the dimension of h. As shown in the appendix 



2 s ( 1 )/ft = 7 /(u;oro) 3 



(30) 



where 

7- 8ae c 3 k ,,2 /(n' 2 h). (31) 

The value of 7 depends on the hber parameters n' 9 and a, cross section of the fiber. Let us consider 10/im area 
silica fibers with n' 2 = 2. 610~ 22 (m/V) 2 as given in (Pil), ([la])- With the values of coefficients given in the appendix 
in MKS units, the coefficient 7 is about 3.610 s . 

With equations P5j) - ([3"Tj) , the transmission coefficient 



T = exp — [- —ln(K c /n)} 7 (32) 

( w Tor 

or 

T=(— )c^F, (33) 

K c 

behaves as a power law, that is smoother than the usual exponential in tunneling amplitudes. The tunneling is 
possible when the exponent in equation (|32p is "not too big". In the semi-classical regime considered above, the 
phase of the wave-function is derived by expansion at lowest order with respect to H. This requires that the exponent 
ln{n c / k) is much larger than unity, then the probability of tunneling is obviously weak. When the 
exponent becomes smaller or of order unity, one is in the "pure quantum limit", and the previous derivation is no 
more valid, since the wave-function cannot reduce to its first order term in H. Nevertheless we can assert by continuity 
argument that tunneling continue to exist, and that it is likely much more efficient. The boundary between these two 
limits can be defined by 



—ln(n c /n) = 1. (34) 

(w To)' 3 

This dependence is drawn in Fig. 5 for the value of 7 given above. The quantum regime is reached as soon as the 
pulse duration is longer than a ps. Therefore quantum tunneling seems within reach of present days experiments. 

To (ft) 
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FIG. 5: Boundary between the quantum and semi-classical regime, for a silica core fiber. The quantum regime stays above the 
frontier. 



C. Quantum switching. 



To estimate the typical length needed for the soliton to tunnel from fiber to fiber, we reason as follows. We estimate 
first the time scale for the quantum tunneling. We split the wavefunction into the 'right' amplitude, and the 
left one, <&£, each one being for the state in one fiber only. Because of the tunneling those states are not eigenstates 
but split into two eigenstates, one even (the ground state) <&s and the other odd, $a, under permutation of the two 
fibers. One has 

$l = (#s ~ $a)/\/2, and <i> R = ($ s + $a)/V2. 

The energy difference between the symmetric and antisymmetric state gives, via the Planck-Einstein relation, the 
typical tunneling time. Let A be half of this energy difference. If at time zero the soliton is on the right fiber, the 
evolution of its amplitude later on is given by 



V 2 

Therefore the amplitude in one fiber oscillates with the period 



(35) 



T osc ~ h/2A. (36) 

In the following derivation, we approximate the energy splitting in each well by using the standard result for a 
particle of momentum p{x) in a double-well: 



2 A = — exp(-- / 

7T « J -a 



\p\ dx) (37) 



where [—a, a] is the x-range under the barrier for the given energy, and u the pulsation of the wave-function in the 
bottom of the well. For a quadratic potential V(x), of curvature V" around the minimum xm, the pulsation of a 
particule of mass m is such that 



V" = moj 2 (38) 

The mass of the particle is deduced from its momentum under the barrier of height Uq, at x — 0, where 
Pq/2to = Uq. Therefore the pulsation writes 



« = y- ir (39) 

where all quantities are in physical units. Note that the dimensions are [P][x] = [S] and [V] = [V]/[x 2 ], therefore 
the dimension of x plays no role. For the fiber problem, we shall consider only one set of conjugate variables, (x, 'J), 
neglecting the rj dependance of the potential, which plays a secondary role, moreover the physical quantities in equation 
(|3"5)) have to be expressed in terms of the reduced ones, and the "time" period of equation (f3"6")l becomes a spatial 
period. This writes 



' Po = S (1) *o 

V" = W^v" = uqS^v" 
< U = W (1) AV 
Tosc — Zk /ro 

where the curvature V" and the height AV of the potential barrier are deduced from equation (fl9|) . that gives 
Al/ = 0.037 + 6k 2 , and v" = -^k~ 2 . Moreover the numerical results in Fig. 4 give = ln(n c / k). In the above 
relations Z is the "true" spatial period along the optical axis of the fibers, obtained from equation (|A1[) after dividing 
all terms by k /ro to obtain a soliton of half-mass equal to unity as assumed in the present section. 
With these expressions, the equation ([36]) becomes 



uj Tosc = Z = cxp S/h, (40) 

to v2w AV 

with S/h= (a J- o) 3 ln(K e //t). 

Since the probability of finding the soliton in a given fiber oscillates with respect to the spatial variable z with a 
wavelength it also oscillates in time from one fiber to the other with the period r = nZ/c, at a given z. Using the 
numerical values given in the appendix for standard fibers, the period of the switching depends on the two parameters 
k and tq. The frequency v = c/nZ, and the spatial period Z are drawn in fig.Q, as function of the coupling 
parameter ratio k c /k. We have chosen two values of pulse duration, To = 0.6ps (dashed line, corresponding to r = 1), 
and l.Sps (solid line, r = 0.1), which are respectively below, and above the frontier drawn in Fig. (5). More precisely 



the dashed line stands into the " semiclassical" regime as soon as k c /k is larger than few units where the WKB 
approximation is valid, whereas the solid line corresponds to the "purely quantum" regime. The two lines displays 
high frequencies, ranging from hundred of Mhz, towards tens of Ghz, that could be interesting for applications to 
high speed transmission. Note that while the solid line corresponds to the pure quantum regime, where the WKB 
approximation used here is not valid, we infer that it could be possible that going beyond the WKB approximation, 
would lead to even higher frequencies. It could then lead to shorter switching lengths than those displayed in Fig. 
d6]-b) ) . In the semi-classical regime, the switching length Z is longer, nevertheless it is much shorter than the half 
period of switching in the CW linear case, L c = 2 ^ hy3 ■ Indeed a pulse duration r = 0.6ps, and a silica fiber, one 
has Le> = 16m, that gives L c = 125— when using k c = 0.2 (see Fig. (4)). For — = 100, the linear half period is 
L c = 12.5km , which is several order of magnitudes longer than the semiclassical switching length Z = 3m (dashed 
curve). 



logio(v) 




FIG. 6: (a) Frequency of the periodic switching v — 1/r, in Log scale, with v in Hz; (b) Spatial period Z, in m, as function of 
the ratio k c /k, for the case to = 0.6ps (r=l, dashed line) and 1.3ps (r=0.1, solid line). For ^ = 100, the period is Z = 10cm 
for the solid line, and about 3m for the dashed line . 



V. SUMMARY AND DISCUSSION 



Even though the tunneling phenomenon is very familiar in many wave-propagation problems, where the "true" 
wave-vector k becomes i k after passing under a classical barrier, (as in the case of evanescent waves in the Fresnel 
theory) , it appears in the present context in a slightly unusual form: starting from the classical model (0 , ([7]) for the 
field enveloppe, which looks strangely similar to the Schrodinger equation, our treatment based on the approximate 
trial functions leads finally to the Euclidean system (|26[) which is not the Schrodinger equation for a wave function. 

Within the trial functions approximation, the WKB or quasiclassical limit gave us the possibility of estimating 
rather easily the rate of quantum tunneling of a single soliton from one fiber to the other, even though it should 
remain classically in the same fiber forever. We found that this rate of tunneling is not small and could well be within 
reach of present day-experiments. 

In the frame of the WKB approximation, we are trying to extend our results by getting rid of the trial functions 
approximation. Our aim is to check if the relation (|29p . that has been shown here to survive when going from two to 
four unknown parameters in the trial functions approximation, is valid beyond this approximation. The calculation 
is heavier than the one presented here, because the time t is now considered as an infinite dimensional parameter, 
then the semi-classical trajectory must be calculated from a set of 4 coupled PDE's, in place of the 4 coupled ODE's 
(|26|) solved here. To derive these PDE's, we can choose for example (^s(u),^s(v)) and (3?(u), $l(v)) as set of conjugate 
variables (p,q) for the classical system ([5]), ©, with Hamiltonian 



n 
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dt 



1 du 



2 1 dt 1 



■ \v\ A — k (uv + vu) 



(41) 



The Euclidean version of equations ©, Q is then obtained by changing (z, $s(u), 3(i>)) into (iz, i$s(u), i$J(u)). The 
important point is that the tunneling factor does not depend on the choice of (p,q), while the Euclidean system 
obviously does. Finally, the heteroclinic Euclidean trajectory is the solution connecting two of the classically permitted 
orbits, from z = — oo to z = +oo. For a given energy E, these are defined by the integro-differential equation 
7~Lnls = E, where ^s(u) = 3(u) = 0, and correspond to the two asymmetrical solitons. 

Because we have found that the Euclidean action S^ hys can be of order of H or even smaller in realistic experimental 
conditions, it could even happen that the WKB quasiclassical approximation is not valid anymore for computing the 



rate of transfer from one fiber to the other. Usually the order of magnitude of the action involved in the soliton 
picture, even in a single fiber, is tacitly assumed to be far bigger than h, which is an assumption distinct from the 
one of a soliton made of many photons. Indeed the soliton picture addresses perturbations to this 'bound state' of 
many photons that may be small enough to imply variations of the action of order of h, and so require some sort 
of ('second') quantization. We plan to come to this general question in future work, and outline here some of the 
estimated problems. 

A treatment using the trial function, but valid beyond the WKB approximation, is obviously more complicated 
than what we did here, and perhaps questionable. Indeed it needs to consider both the trial functions and their 
parameters as operators. Moreover it amounts to assume that the fluctuations in t and z are decoupled, and, last 
but not least, our result derived in the WKB approximation likely signals that the assumption behind the classical 
(meaning non quantum) theory for describing soliton in coupled fibers does not hold anymore and that the quantum 
picture has to be used from the start, which makes it theoretically challenging. 

We assumed that every phenomenon under study involved solitons seen as a coherent quantum objects. We argued 
that this requires that any typical time, the tunneling time in particular, is far shorter than the coherence time. 
This coherence time is of order of t c /N , with t c coherence time of a single photon in the soliton, i.e. its mean-free 
flight time without change in phase or frequency. Because of the division by N this may be a very short time. At 
times longer than the coherence time any physical effect related to the quantum coherence between states of solitons 
propagating in either fiber is washed out. The final state, as described in the density matrix formalism, is a state 
of equal probability of the soliton on either side without nondiagonal element. The experimental manifestation of 
this state will be a probability 1/2 of observing a soliton in cither fiber without any possible interference between 
the states on either side. Somehow this will bring the system back to a fully classical state, except that this classical 
state has a probabilistic underpinning that is absent from the classical system: in the fully classical system the soliton 
remains always in the same fiber, although in the quantum one its final state has a probabilistic nature. 

Looking at the other side of the coin one realizes that, because the soliton is a composite object, and if it remains 
coherent during a sufficiently long amount of time, its phase is the phase of a single photon multiplied by the number 
of photons. Therefore any interference experiment between coherent soliton states will have much narrower interfringe 
than with a single photon or incoherent photons, this interfringe being the one for a single photon divided by the 
number of photons making the coherent soliton. This could be of interest for gyroscopes based on the Sagnac effect 

m- 

Indeed a central issue concerning the observability of the tunneling effect we present in this communication is the 
one of the quantum coherence of the soliton, related itself to all dissipative effects that can break up this coherence, 
and that makes the main topic discussed in the present special issue. Nevertheless, even if the coherence is limited, 
there is still tunneling, but at a reduced rate In that case we suggest to use twin fibers with coupling coefficient 
k periodically modulated in z, in order to stimulate the switching process. 
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APPENDIX A: PHYSICAL UNITS 



This appendix is about the relationship between quantities measured in physical units for a standard fiber carrying 
solitons and the dimensionless quantities used in the bulk of our paper. We relate first the number of photons in a 
typical soliton to its time duration, a duration called tq that we shall use afterwards to give various order of magnitudes 
pertinent to our problem. 

We use the standard expression of the electric field in a fiber, written with the same notations as in the book by 
Newell and Moloney ■ The electric field of the EM wave in the fiber is modelized by the wave-packet expression, 
E = R(xt)B(z, t) exp(icjoi — koz) + c.c. which obeys in a first approximation Maxwell's equations, when the duration 
of the pulse is not too short, x± being the transverse coordinate, and R(xt) the dimensionless radial amplitude with 
J \R(x')\ dx' — (7 as the core area. 

Taking ^ as a small parameter, Sto being the frequency width of the pulse, one obtains the NLS equation written 
in variables z, r = t — z/v g , the nonlinear and the dispersion term having opposite signs: 

dB kl d 2 B 2Trn' 9 , l2 , A , 



Note that we turned to the standard writing of the coefficients of the NLS equation, n' 2 being the modulus of the 
coefficient of the cubic Kerr effect and fc the modulus of the second derivative of the wavenumber fco with respect to 
the the frequency of the EM wave. 
The soliton solution is 

B(z, t - z/v g ) = 5 m sech ((t - z/v g )/r ) , (A2) 

with 

f3\B m \ 2 = l/rl (A3) 



where f3 = Its energy is 



W = J P(t)dt = a J I(t)dt, (A4) 

where I(t) is the optical intensity measured in watt per square meter, and P(t) is the Poynting vector integrated 
across the fiber section, with a result expressed in Watts: 

P{t) = J (EAH).zdS (A5) 

where z is the unit vector in the direction of propagation. 

The magnetic field in the wave (supposing it is linearly polarized with the electric field in the x-direction) is 
H x = nceoEy, n index of refraction. For a material with instantaneous response, after averaging over one period of 
the field oscillations, one finds 

P(t) = 2ce an\B{t)\ 2 . (A6) 
At leading order, i.e. by taking into account the linear part of refractive index, n = hq, this yields 

I^(t) = 2cs n \B(t)\ 2 . (A7) 

Whence the energy of the pulse is: 



f A k 1 

W (0) = Nhuj = 2n c£ cr / \B(t - z/v g )\ 2 dt = 4n ce <T-^- — . (A8) 
J 2-Kn 2 t 

Finally the relationship between the photon number and the pulse duration writes 

NTQ = Ano£oa ^\ (A9) 

Similarly the action is 

5 (0) =Nh = W (0) /lu, (A10) 

at leading order. 

In MKSA units, with standard values (see [ll|) of optical fibers composed of silica cores, this gives: 
no = 1.5 

eo = 0. 89. 10 _ 11 F/m, or cqc — 1/Zq with Zq = 3770 the impedance of free space, 
h = lCT 34 J.s 
A = 1.55.10 _6 m 

a ~ 10 -11 for a 10/im 2 -area fiber, 
fco = 2.2.10- 26 s 2 /m 
n' 2 = 2.6.10" 22 (m/V") 2 

With these data, the number of photons in the pulse of duration tq, measured in seconds obeys the relation: 



Nt ~ 3.1(T 5 , 



(All) 



that gives N ~ 3.10 7 photons for a ps-pulse. 

Let us note that the nonlinear index of refraction n' 2 may be several orders of magnitude larger, when using other 
materials. For example, in the experiment of Wa et al. [ill ], the optical switch was studied in multiple quantum well 
wave-guides, with n 2 = 10~ 13 (m/V) 2 . 



Coherent part of the energy and action 



Let us write the energy and action as 

W = W {a) +W (1 \ (A12) 

and 

S = S {0) + S {1) . (A13) 

The dominant contributions are proportional to the linear part of the refractive index, namely a term contained 
in the Maxwell equation. The subdominant contributions and , correspond to the terms contained in the 

envelope equation, they are perturbations to the dominant effects calculated above. For two coupled fibers, these 
perturbations result from balanced effects of dispersion, coupling and nonlincarity. They are proportional to the 
energy and action of the the dimensionless NLS equation ((6][7j) , 

I - ^ (Ai4 » 

where the scaled energy 7Ynls is defined in equation (|41[1 , and the action in equation pop. 

The coefficient u>W may be calculated by using the expression of the Poynting vector (|A6p valid for dispersionless 
Kerr media, where 

n = n a + n 2 I = n + n' 2 \B(t)\ 2 , (A15) 
This gives W^' Kerr ^ = 2n' 2 ce a f dt\B(t)\ 4 . 

Taking the hyperbolic secant solution (|A2IIA3[) , one obtains W^ 1,Kerr ^ = 2n' 2 ceo<7 \B M (t)\ 4 T Q J dtsech\t), or 



W iX,Kerr) = .4^^ W ^S, (A16) 

P T 

where Hnl" i s ^he Kerr contribution of the Hamiltonian (second term in the r.h.s. of equation [5|) . This correction is 
the " coherent" part of the energy in the sense that it is proportional to the square of the intensity, or of the photon 
number. Finally we are ready to express the physical value of the Hamiltonian and action associated to a soliton 
whose temporal width is scaled to tq as it was assumed in sections 3-4, by using the expression (fAT4| with 

w& si 4n 2 ce a-^ . (A17) 
P T o 

This allows in particular to express concretely the constraint that nonlinear effects are small, that is that << 
W(°\ a condition equivalent to 

2n' 1 , A nN 

^ = 3& <<L (A18) 



Finally the physical action associated to the quasi-classical trajectory is approximately given by the expression 

s phys = ^_g E ( Al9) 

The quantum tunneling coefficient T — exp(—2Sg lys /h) is expected to be experimentally observable when the 
physical action is not too large with respect to h. When Sg ya becomes of order h, we can infer that we are in the 
"pure quantum" regime. The frontier between these two regimes may be drawn in the space parameters {k,tq) by 
introducing the parameter 



that allows to write 



7 . 8 y° (A20) 



2 «^) = 7 (A21) 



With the data given above for standard fibers, 7 = 3.610 8 , the semi-classical regime stands below the curve drawn 
in section HVl (see Fig. (5)). Consequently pulses longer than one ps typically stands in the pure quantum regime. In 
conclusion we predict quantum tunneling for realistic conditions of soliton propagation in two coupled fibers. 
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